Free flux flow magneto-resistance of single crystal BSCCO 
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Measurement of the free flux flow resistance in monocrystalline highly anisotropic BSCCO reveals 
a low field magneto-resistive effect nearly 100 times that expected from a naive application of the 
Bardeen-Stephen rule R/R„ — B/B C 2- Furthermore in the vortex solid phase it saturates to a 
constant value at about i/ C 2/100 as if the vortices were moving 100 times faster than expected. An 
attempt to account for this behavior by a combination of sheared planar vortex flow obeying the 
BS rule and the quasi-particle c axis conductivity between the Josephson coupled superconducting 
planes is only partially successful. In particular it is unable to account for the saturation feature 
which occurs in both optimally and underdoped samples in the low temperature vortex solid phase. 

PACS numbers: 74.72.Hs, 74.25. Fy, 74.25. Qt, 74.25.Sv 
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I. INTRODUCTION 

High T c cuprate superconductors differ from their low 
T c counterparts in their extreme type II character with 
a Ginzburg-Landau parameter k ~ 200 corresponding 
to extremely short coherence length £ ~ 1 nm and high 
upper critical fields B C 2 ~ 100 T. The small Abrikosov 
vortex core size of order £ results in a series of bound 
states with energies situated in the gap whose mean en- 
ergy spacing is of order < Se >~ A 2 /Ep ~ h 2 /mt; 2 
where A is the superconductor gap energy (modulus of 
the order parameter), Ep the Fermi energy and m the 
effective mass. The core structure can affect the mobility 
for translational motion of the vortices by modifying the 
spectral properties of the electron momentum transfer. 

Dissipation in type II superconductors in magnetic 
field H > H c i in slowly varying electromagnetic fields 
is dominated by the dynamics of Abrikosov vortices [l[. 
Material details of the superconductor enter mainly via 
the "friction coefficient" that gives the vortex velocity in 
terms of the driving Magnus-Lorentz force. The friction 
is controlled by the spectral density and relaxation prop- 
erties of the low-lying "core states", i.e., quasiparticle 
states localized to the vortex core. In dirty supercon- 
ductors, where the mean free path £ <C £, the friction 
coefficient is related to the normal state resistivity by 
the Bardeen-Stephen (BS) rule 0: 



PS = ap n B/B c2 



(1) 



where p/ is the flux flow resistivity arising from vortex 
motion in the absence of pinning, p n is the normal state 
resistivity, B is the magnetic field, B C 2 the upper critical 
field and a — 1. A more careful examination of the core 
states allowed the extension of Eq. |1| with a « 1 to 
the moderately clean limit £ > £ as well [l|. The BS law 
has been experimentally confirmed for a broad range of 
conventional (s-wave) superconductors Q. 

In unconventional superconductors with gap nodes, 
BCS theory suggests a high density of core states at the 



Fermi energy. Nevertheless calculation [J] of the flux- 
flow resistivity revealed that as long as one stays in the 
moderately clean limit, Eq. ([T]) remains valid with a 
possible enhancement of a but still of order 1. Recent 
measurements Q on several anisotropic non-high- T c su- 
perconductors that most likely exhibit gap nodes confirm 
Eq. ((T|) with a moderately enhanced 1.6 < a < 4.7. 

The problem of low-temperature vortex dynamics in 
BSCCO (Bi 2 Sr2CaCu208+(5) is particularly interesting 
in this respect in the light of the recent results on the 
structure of the vortices in this material. Scanning tun- 
neling microscope (STM) spectroscopy @ revealed that 
the zero-energy peak in the density of core states is miss- 
ing. This, together with the results of inelastic neutron 
scattering Q and NMR § experiments, led to the sug- 
gestion that some concurrent, non superconducting, or- 
der exists at the vortex cores when the superconducting 
order parameter is suppressed Q. The idea received di- 
rect support from STM spectroscopy 10] where a peri- 
odic modulation of the local density of electronic states 
around the vortex cores was observed. Since the decay 
length of this modulation is much longer than the su- 
perconducting coherence length, not only the structure 
of the core, but also the structure of the flow field for 
vortex transport should be different from conventional 
Abrikosov vortices with probable consequences for the 
velocity-force relation. 

The experimental situation in the high-T c supercon- 
ductors, and in particular in BSCCO, is not clear. Low- 
frequency transport measurements [111 ] in the vortex liq- 
uid phase close to the critical temperature T c are in rea- 
sonable agreement with the BS law. Microwave and mil- 
limeter wave impedance measurements [l2j at low tem- 
peratures, on the other hand, indicate a large and nearly 
£?-indcpendcnt dissipation in a broad field range. 

To clarify the situation and to see if these very spe- 
cial aspects of the vortex structure influence the dynam- 
ics, we have performed a systematic investigation of the 
flux-flow resistance in BSCCO single crystals by measur- 
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ing the a-b plane voltage-current {V-I) characteristics up 
to currents well above the threshold current for dissipa- 
tion in c-directed magnetic field. The free flux- flow re- 
sistance as measured on an ab face can be approximately 
described as oc B 1 / 2 for low fields followed by saturation 
(becomes field independent) above about f T in the low- 
temperature vortex solid phase. As the saturated value 
corresponds quite well with the extrapolation of the nor- 
mal resistance to low temperatures, it might be naively 
interpreted as reflecting a ~ 100 in Eq. (fTJ), at least at 
1 T, indicating that vortices move at up to two orders of 
magnitude faster than predicted by BS law in low fields. 
Although it is known that vortex velocities much larger 
than the BS value may result from a nonlinear instability 
at high vortex velocities [HI], not only does its onset in 
BSCCO films[ll| occur at current densities about two or- 
ders of magnitude higher than investigated in this study, 
the characteristic signature (TH . Ibj ] of a nonlinear run- 
away in the V-I curve is absent in our case. 



II. EXPERIMENT 

Measurements were made on nine single crystals from 
three different batches of Bi2Sr2CaCu20s+,5 fabricated 
by a melt cooling technique [la]. The typical dimen- 
sions of the crystals were 1 x 0.5 x 0.003 mm 3 with the 
shortest dimension corresponding to the poorly conduct- 
ing c axis. Most of the crystals have been close to op- 
timal doping with a resistance-determined critical tem- 
perature T c « 89 K and transition width about 2 K in 
zero field. The diamagnetism in a f Oe field set in in a 
uniformly progressive way below T c to near 100% at low 
temperature. Two underdoped crystals exhibited a resis- 
tive T c w 51 K and T c 53 K with width 6-8 K and slow 
initial diamagnetism onset to about 30 K before proceed- 
ing to about 40% flux exclusion at low temperature. 

The resistance measurements were made in the usual 
four point configuuration on an ab face of the crystals. 
The contacts to inject and withdraw the current are ei- 
ther on the top a-b face or encompass the ends of the 
crystal extending to both opposing a-b faces. The fea- 
tures reported here are common to both geometries. Four 
more contacts on the same a-b face as the current injec- 
tion serve to measure the voltage parallel and transverse 
to the current flow. The contacts are made by bonding 
25 fim gold wires with silver epoxy fired at 900 K in an 
oxygen atmosphere resulting in contact resistances of less 
than 3 f2 for the current contacts. The sample is inserted 
in the bore of a superconducting magnet with the field 
along the c direction. 

To measure at the high currents required for the free 
flux flow regime without significant Joule heating, we ap- 
ply short (typically 50 /xs or less) current pulses of isosce- 
les triangular shape with a repetition period of 0.2 to 1 
s. Further technical details are given and the issue of 
Joule heating for the such experimental conditions is an- 
alyzed in Ref. [161 ] with the conclusion that the temper- 
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FIG. 1: Typical voltage-current characteristics in B — 3 T 
at temperatures (from left to right) 60, 45, 38, 19, 11, 7, and 
5 K. Inset: Magnified view of the low-current region. The 
dashed lines are linear fits. (See text for the definition of I a f, 
Ith, Rf, and R th .) 



ature change in the area between the voltage contacts is 
negligible for the duration of the pulse. 

Typical V-I characteristics measured at different tem- 
peratures and in a field of 3 T arc shown in Fig. [T] At 
temperatures below the vanishing of zero current resis- 
tance, which we interpret as the freezing of the vortex 
system (in B = 3 T this happens at T m — 33 K), dissipa- 
tion sets in abruptly beyond a sharply defined threshold 
current I t h- Usually, but not always [l6| . a linear seg- 
ment appears above the threshold current (see insert) 
the slope of which we denote by R t h- At higher currents, 
the differential resistance increases further and saturates 
at a value Rf S> Rth in the high-current limit. The high- 
current linear segment of the V-I characteristic extrap- 
olates to finite current If ^> Ith at zero voltage. Above 
the melting temperature, a finite differential resistance 
Ri^o is observed in the I — ► limit. The V-I curve, 
however, is still nonlinear and I s f is finite up to a tem- 
perature T* situated between the melting temperature 
and the critical temperature T c . 

The temperature dependence of the low-current and 
high-current differential resistances, Ri->o and i?/, are 
shown in Fig. [2]for two magnetic fields. Between T c and a 
field dependent characteristic temperature T* , these two 
differential resistances are identical (the V-I curve is lin- 
ear). At T* the two resistance curves bifurcate; Ri^o de- 
creases rapidly with decreasing temperature and reaches 
zero (< 100 fiCl) at what is taken to be the melting tem- 
perature, T m , while the high-current resistance Rf re- 
mains finite and varies smoothly across the melting line. 
The "pinned liquid" domain extends over T m < T < T* 
where signs of pinning are still present. Below T m , in 
the vortex solid phase, Ri^o is zero and Rf attains a 
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FIG. 2: Temperature dependence of the low-current (fl/_»o) 
and high-current (Rf) differential resistances in 1 and 3 T 
magnetic fields. Ri^o is also shown for zero field for reference. 
(See text for the definition of T m , T c , and T* .) 
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FIG. 3: Magnetic field dependence of the high-current Rf 
differential resistance at different temperatures. Insert shows 
the low-current (Rth) differential resistance. 



finite, temperature-independent value visible in Fig. [2] 
Moreover, this value is the same for 1 T and 3 T and 
agrees reasonably with a linear extrapolation of the nor- 
mal phase resistance to low temperature. 

More insight on the low-temperature saturation of the 
high-current resistance Rf is gained from its field depen- 
dence for several temperatures shown in Fig. [5] In the 
liquid phase, we find an approximately B 1 / 2 field depen- 
dence over the whole field range. In the solid phase, how- 
ever, Rf becomes independent of magnetic field above 
about 1 T. On one sample the measurements were ex- 
tended to 17 T at 5 K and no variation of Rf in excess of 
the experimental uncertainties of about 10 % were found 
in the field range 1 to 17 T, in contrast to the 17-fold 
increase predicted for the ab resistivity by Eq. fTJ). The 



differential resistance close to the threshold current, Rth, 
behaves similarly (see inset). In fact, in the parameter 
range where both quantities can be measured, we find 
a field and temperature independent proportionality be- 
tween R t h and Rf. 

The behavior of the characteristic currents 7th and If 
further corroborates the similarities of the dissipation 
mechanisms close to threshold and in the high-current 
limit. The similarity is again underlined by the same 
B~ x / 2 field dependence for both quantities in the low 
temperature region. All these observations suggest that 
the same mechanism is responsible for the dissipation 
in the vicinity of the threshold current as in the high- 
current limit. The similarity of behavior is attributed 
to the way in which the resistive front propagates from 
the current contacts towards the middle of the sample 
as the current is increased 17J • The current distribution 
reaches the resistive limit characterised by Rf only after 
the two resistive fronts meet. As described earlier[l6j at 
low temperature there is a difference between field cooled 
(FC) and zero field cooled (ZFC) preparation. The ZFC 
data for both I t h and If exhibit a characteristic peak at 
the same line in the (B, T) plane which defines a low 
temperature region where the FC prepared state is only 
metastable. Both FC and ZFC preparations give rise to 
the same Rf. 

A first analysis to get our bearings consists of nor- 
malizing the resistance to the resistance measured in the 
normal phase and the field to B c2 . A refinement on this 
is to estimate what the normal resistance would have 
been at the temperature of the measurement were the 
sample not superconducting. Following this procedure 
and noting that when the magnetoresistance value satu- 
rates with field its value approximates the expected nor- 
mal resistance, it has been hypothesised to be just that 
and the value used to better interpolate the normal re- 
sistance. For the upper critical field, we use the form 
B c2 (T) = (120 T)[l - (T/T c ) 2 } yielding dB c2 /dT\ Tc = 
-2.7 T/K pi]. To estimate R n in the T < T c range, we 
use a quadratic interpolating function to connect with 
the T > T c range on the hypothesis that the field satu- 
ration value for Rf corresponds to R n {T). For the very 
low field data below about 0.5 T where we could not ap- 
ply sufficient current to reach the Rf asymptote, we used 
the temperature and field independent scaling factor be- 
tween Rf and Rth found at higher fields to estimate Rf 
from the measurement of Rth- The flux flow resistance 
data treated in this way are displayed in Fig. 0J 

It is evident from Fig. 2]that the ratio of flux flow resis- 
tance to normal resistance is very different from the field 
proportional BS p a b resistivity ratio of Eq. ([!}, ranging 
from a y/~B like variation at low field to a saturation value 
for the low temperature vortex solid phase. Nonetheless 
there is a trend that with increasing temperature, in the 
vortex liquid phase, the saturation feature disappears, at 
least for the field values accessible to us. In the solid 
phase at low temperatures there is strictly no field de- 
pendence of Rf whatsoever above about 1 T although as 
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FIG. 4: Reduced differential resistance (Rf / R n ) as a func- 
tion of the reduced magnetic field (B/B^). The BS law for 
the flux flow resistivity is represented by a straight line and 
corresponds to the right hand scale. 



mentioned above its value is in good agreement with the 
linear extrapolation to T = of the normal resistance 
measured above T c . Since the field where saturation oc- 
curs is in the order of 10~ 2 B C 2, it might appear that for 
fields below saturation, the vortices move ~ 10 2 times 
faster than expected, or that the vortex friction coeffi- 
cient is ~ 10 2 smaller, from a simple application of the 
BS relation. 



III. FRAMEWORK FOR UNDERSTANDING 

The experiment measures resistance and the ratio of 
resistances is only proportional to the ratio of resistivities 
if the current distributions are identical. This is not the 
case in anisotropic materials like BSCCO where the cur- 
rent penetration depth is determined by the anisotropy of 
the resistivity. The latter is not only different in the su- 
perconducting state but is also field dependent. From 
this point of view we might interpret the comparison 
made in Fig. 3] as an indication that the current pene- 
tration is considerably less in the low temperature vortex 
solid phase, progressively approaching the normal phase 
penetration depth at higher temperatures in the vortex 
liquid phase. Indeed in the liquid phase, the data are in 
rough agreement with BS law if the effects of inhomoge- 
neous current distribution within the sample are taken 
into account. 

The voltage measured on the top plane of the sam- 
ple is affected by the c-axis properties because they in- 
fluence the distribution of the transport current within 
the sample P3 ■ In the simplest approach we could sup- 
pose that well into the ohmic response regime the cur- 
rent is distributed as for a normal anisotropic ohmic 
conductor [l9j. a surmise which is borne out by numer- 



FIG. 5: Schematic illustration of voltage-current response 
between Josephson coupled planes. Part (a) shows a stan- 
dard SIS junction between s-wave superconductor planes, 
(b) shows d-wave response as found in the experiment of 
Ref. [2(il] where the d-node quasi particle shunt conductance 
a c qv m cr^/40. (c) shows the effect of spatial dephasing across 
the junction. 



ical simulations based on the superconductor model of 
Ref. [13] • The penetration depth d for a length I between 
current injection points on the surface is then given by 
d 2 w (a c /a ab )£ 2 (provided d <C sample thickness) and the 
resistance R along the length of the sample will scale like 
- ab (£/d)w 



R « p ab (£/d)w w wy/p ab p c . If p ab = pf of Eq. and 
the resistive response of the Josephson coupled planes in 
the superconducting phase is represented by p c s , one ob- 
tains R f /R n = (p c J 'p^/Zy/aB/Ba. Results close to T c 
are indeed well described by this form. Putting p c s = p c nl 
the 75-K data in Fig. [4j for instance, are reasonably 
well fitted with a = 4, which, given the uncertainties of 
the simple resistive thick-sample estimate, we regard as 
rough agreement with BS law. The crux of the matter 
then, from this viewpoint, is in p^(B,T). 

The c axis conductivity a c s can be obtained from the 
I — V characteristics of the interplane Josephson junc- 
tions. The usual representation of a Josephson junction 
in a classical s-wave superconductor subject to a spa- 
tially uniform phase difference is illustrated on Fig. [5^l: 
the voltage response to current is nil up to the Josephson 
current Jj at which point the junction opens to V = V„ = 
2A/e and responds to further current according to the 
normal phase tunnel conductance a n = 1/ p™ which is re- 
lated to the critical current by the Ambegaokar-Baratoff 
relation Jj = ircr n A/2e — (ir/4)a n V g . For currents 
J 3> Jj the c axis resistivity is that apppropriate to the 
normal state p% — p c n . Assuming the penetration depth 
to be limited by the resistive anisotropy gives the result 
quoted above, but this regime would only be attained 
for extreme currents ~ 10000 Acm~ 2 corresponding to 
I ^ 20 A for our samples. That is clearly not the regime 
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which concerns us here. 

The answer to this problem seems to be that the 
Josephson critical currents are much lower than the 
Ambegoakar-Baratoff prediction. Measurements [2Cj on 
micrometer size mesa stacks of BSCCO junctions at zero 
field show much lower critical currents, Jf ~ (n /2)<j qp V g 
where the quasi-particle conductivity a qp ~ cr n /A0 as 
V — > at low temperature and V g ~ 50 mV yield- 
ing J j w 500 Acm~ 2 . As illustrated schematically on 
Fig. \5\3 the opening of the junction is followed by a 
voltage response at higher currents with slope p n when 
J > J j 3> J j- But even this intermediate regime is not 
observed: a penetration depth of 150 nm correspond- 
ing to an anisotropy factor 7 ~ 3000 would involve 100 
planes and a voltage response of 5 V, still a factor of 
~ 10 2 higher than observed in our measurements. To 
understand what we see, we must take account of the 
non-uniformity of the phase across the planes due to 
the random positioning of vortices from plane to plane. 
The critical current for the junction is then reduced to 
JJ = (cosa(r)) jj where a(r) represents the phase dif- 
ference across the junction at position r in the plane and 
(cosa(r)) the value of cos a averaged over the junction 
area on a Josephson penetration length scale [U [22J. It 
is typically of order 10~ 2 . As illustrated on Fig. 
the junction will then open to a potential difference of 
(cosa(r)}V g <C V g due to phase slippage at J = JJ, 
beyond which it should continue to see a tunneling con- 
ductance a qp of quasi-particles in the gap until V > V g , 
when the tunneling of quasi particles at the gap edge will 
result in the usual normal phase dynamic conductance. 
The gap quasi-particle conductance is thus experienced 
over a wide current range between JJ and JJ j (cosa(r)) 
over which the resistance ratio is considerably enhanced 
to Rf/Rn = (<jfJ<7 c qp ) 1/2 VaB/B c2 where (a c Ja c qp ) « 
30 — 40 at low temperature [20(, again if we assume Eq. 
([I]) to be valid for the second square root term. The 
magnetic field dependence of a qp ~ cr qp (l + 0T)B/B c2 ) 
was measured in a subsequent experiment [23j, also on 
micrometer size mesa samples. If finally we include also 
the effect of quasi-particle conductance a qp (B) in the ab 
plane, we find 

Ik = [K / aB/B c2 

Rn V °w V (1 + W^X 1 + (af p /a-")B/B c2 ) [ 1 

This is the expression that we compare with the ex- 
perimental results on Figs. E^a)-(d). 

IV. INTERPRETATION OF RESULTS 

Equation {2]| is compared with the experimental find- 
ings for both the optimal doped and the underdoped sam- 
ples on Figs. [6] For simplicity a single plot of Eq. ([2]) 
is drawn as a dashed line on the log-log plot against 
B/B c2 (T) where B c2 (T) is the temperature corrected 
value as described above. It corresponds to a choice of 



parameter values appropriate to the optimally doped sit- 
uation at low temperature T -> 0: B c2 (T = 0) = 120 T, 
a c Ja c qp = 300, a = 1,0= 10, af p /a a n b -> and 
R n /R n (T+) = 1/3 corresponding to a linear extrapola- 
tion to T = . A parallel solid line is also drawn as the 
best \/B/B c2 fit to the data. 

We note that: 

1. The order of magnitude is correct to within a factor 
of 2 for both optimal and underdoped samples, with the 
notable exception of the saturation plateau values in the 
vortex solid phase, especially in the underdoped sample. 

2. The square root behavior predicted by the BS re- 
lation is approximately correct, again with the obvious 
exception of the plateaus. 

3. The magneto- resistance plateaus are all situated in 
the vortex solid phase. It is to be remarked that the 20 K 
magneto-resistance in the underdoped sample begins to 
saturate, but on reaching about 0.5 T it reverts to an 
increasing power law. The melting field for 20 K is es- 
timated to be about 0.3 T from the appearance of finite 
Ri^o- 

Clearly Equation {2} does not describe all the features 
of the magneto-resistance, even if it does reasonably well 
for temperatures above the vortex solid melting and in 
the solid region for fields below the saturation. It has a 
structure which could in principle describe a plateau if 
the (3 coefficient in the field dependence of the shunt con- 
ductance of the Josephson coupled layers were to be much 
larger (~ 10 2 ) such as might be introduced by aligned 
vortex core NIN tunneling junctions or possibly e?-wave 
node quasi-particle to vortex core tunneling. Such ef- 
fects have not however been seen in the mesa magneto- 
conductivity measurements [2^], although these were per- 
formed at relatively high temperature (> 20 K) and high 
field, so in the vortex liquid phase. 



V. CONCLUSION 

It is difficult to conclude that the BS relation is untrue, 
either by a multiplicative factor or in its form. To deter- 
mine the multiplicative factor requires a good way of es- 
timating the normal resistivities in the superconducting 
phase. Also, because the agreement is reasonably good 
at high temperature and the vortex entity is not expected 
to change nature from high to low temperature, the seri- 
ous disagreement at low temperature is much more likely 
to be related to the nature of the solid phase as the BS 
relation is a single vortex property. 

It can however be concluded from the sublinear, ap- 
proximately BJl 2 , field dependence that if the linear field 
dependence of the BS relation is correct, the current dis- 
tribution is not constant and its penetration depth is con- 
trolled by the resistive anisotropy and thus by pj. 



6 




-i 



H H [Tesla] H / H 2 (T) 

(c) (d) 

FIG. 6: Comparison of experimental measurements with Eq. (|2]) for the flux flow resistance as determined by anisotropy using 
the Bardeen-Stephen flux flow relation. The raw experimental results are given in a linear representation (a) and (c) and in 
a reduced field representation which takes account of the temperature variation of B C 2 = B C 2(T) — 120(1 — (T/T c ) 2 T in a 
log- log representation in (b) and (d). The dashed lines represent the formula of Eq. using the parameters indicated in the 
text. The same parameters are used for both optimal and underdoped samples for want of data specific to the underdoped 
case. The solid lines are simply an indication of a (B / B^) 1 ^ 2 power dependence. 
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